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Abstract. In this paper we present algorithms that compute cer- 
tain local cohomology modules associated to a ring of polynomials 
containing the rational numbers. In particular we are able to com- 
pute the local cohomological dimension of algebraic varieties in 
characteristic zero. Our approach is based on the theory of D- 
modules. 



1. Introduction 

1.1. Let R be a commutative Noetherian ring, I an ideal in R and 
M an i?-module. The i-th local cohomology functor with respect to / 
is the i-th right derived functor of the functor Hj(-) which sends M 
to the J-torsion UfeLi(0 '-m I k ) of M and is denoted by H}(—). Local 
cohomology was introduced by Grothendieck as an algebraic analog of 
(classical) relative cohomology A brief introduction to local cohomol- 
ogy may be found in appendix 4 of ||. 

The cohomological dimension of / in R, denoted by cd(R,I), is the 
smallest integer c such that the local cohomology modules Hj(M) = 
for all q > c and all i?-modules M. If R is the coordinate ring of an 
affine variety X and / C R is the defining ideal of the Zariski closed 
subset VCX then the local cohomological dimension of V in X is 
defined as cd(R,I). It is not hard to show that if X is smooth, then 
the integer dim(X) — cd(R, I) depends only on V but neither on X nor 
on the embedding V X. 

1.2. Knowledge of local cohomology modules provides interesting in- 
formation, illustrated by the following three situations. Let J C R and 
c = cd(R, I). Then I cannot be generated by fewer than c elements. In 
fact, no ideal J with the same radical as / will be generated by fewer 
than c elements. 

Let H l dR stand for the i-th de Rham cohomology group. A second 
application is a family of results commonly known as Barth theorems 
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which are a generalization of the classical Lefschetz theorem that states 
that if Y C is a hypersurface then H dR (¥^) — > H dR (Y) is an iso- 
morphism for z < dim(y) — 1 and injective for % = dim(F). For ex- 
ample, let Y C Pg be a closed subset and I C R = C[xq, . . . ,x n ] the 
defining ideal of Y. Then if^(P^) — > H dR (Y) is an isomorphism for 
i < depth (Cy) — cd(R, I) (compare (TTS|, 4.7 and ||, the theorem 

after 111.7. 6). 

Finally, it is also a consequence of the work of Ogus and Hartshorne 
([0, 2.2, 2.3 and §, IV.3.1) that if / C R = C[a? , . . . ,ac n ] is the 
defining ideal of a complex smooth variety V C Pg then, for z < 
n — codim(V), dime socR(H^(H n ~ l (R))) equals dime (V", C) where 
Hl(V,C) stands for the z-th singular cohomology group of the affine 
cone V over V with support in the vertex x of V and with coefficients 
in C (soc^(M) denotes the socle (0 \m m) C M for any i?-module M). 



1.3. The cohomological dimension has been studied by many authors, 
for example R. Hartshorne (0), A. Ogus (| ]15||) , R. Hartshorne and 
R. Speiser (0), C. Peskine and L. Szpiro (fj), G. Faltings (§), 
C. Huneke and G. Lyubeznik (||). Yet despite this extensive effort, 
the problem of finding an algorithm for the computation of cohomo- 
logical dimension remained open. For the determination of cd(i?, /) 
it is in fact enough to find an algorithm to decide whether or not the 
local cohomology module H^R) = for given i, R, I. This is because 
Hj(R) = for all q > c implies cd(R, I) < c (see 0, section 1). 
In G. Lyubeznik gave an algorithm for deciding whether or not 
H}(R) = for all / C R = K[xi, . . . , x n ] where K is a field of positive 
characteristic. One of the main purposes of this work is to produce 
such an algorithm in the case where K is a field containing the rational 
numbers and R = K[xi, . . . , x n ]. 

Since in such a situation the local cohomology modules H](R) have 
a natural structure of finitely generated left D(R, i^)-modules ([|Tl"|j), 
D(R, K) being the ring of .fT-linear differential operators of R, explicit 
computations may be performed. Using this finiteness we employ the 
theory of Grobner bases to develop algorithms that give a representa- 
tion of Hj(R) and H l m (Hj(R)) for all triples i,j E N,I C R in terms 
of generators and relations over D(R,K) (where m = (xi, . . . ,x n )). 
This also leads to an algorithm for the computation of the invariants 
Xi ) j(R/I) = dimft- socr^H^H^ 3 (R))) introduced in fTT| . 



We remark that if R is an arbitrary finitely generated i^-algebra and 
I is an ideal in R then, if R is regular, our algorithms can be used to 
determine cd(-R, /) for all ideals I of R, but if R is not regular, then 
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the problem of algorithmic determination of cd(R, I) remains open (see 



also the comments in subsection |6.4|) 



1.4. The outline of the paper is as follows. The next section is devoted 
to a short survey of results on local cohomology and D-modules as they 
apply to our work, as well as their interrelationship. 

In section |3| we review the theory of Grobner bases as it applies to 
A n and modules over the Weyl algebra. Most of that section should 
be standard and readers interested in proofs and more details are en- 
couraged to look at the book by D. Eisenbud (|§, chapter 15 for the 
commutative case) or the fundamental article || (for the more general 
situation. 

In section |4] we generalize some results due to B. Malgrange and 
M. Kashiwara on .D-modules and their localizations. The purpose of 
sections [| and |5] is to find a representation of Rf <S> N as a cyclic A n - 
module if N is a given holonomic -D-module (for a definition and some 
properties of holonomic modules, see subsection |2.3| below). Many of 
the essential ideas in section | come from T. Oaku's work [fPfl . 

In section ^ we describe our main results, namely algorithms that for 
arbitrary i,j,k,I determine the structure of Hj(R), H^ a (Hj(R)) and 
find Xij(R/I). Some of these algorithms have been implemented in 
the programming language C and the theory is illustrated with exam- 
ples. The final section is devoted to comments on implementations, 
effectivity and examples. 

It is a pleasure to thank my advisor Gennady Lyubeznik for suggest- 
ing the problem of algorithmic computation of cohomological dimension 
to me and pointing out that the theory of .D-modules might be useful 
for its solution. 



2. Preliminaries 

2.1. Notation. Throughout we shall use the following notation: K 
will denote a field of characteristic zero, R = K[x\, . . . ,x n ] the ring 
of polynomials over K in n variables, A n = K(x±,di, . . . ,x n ,d n ) the 
Weyl algebra over K in n variables, or, equivalently, the ring of K- 
linear differential operators on R, m will stand for the maximal ideal 
(xi, . . . , x n ) of R, A will denote the maximal left ideal (dx, . . . , d n ) of 
A n and / will stand for the ideal {fx, ■ ■ ■ , f r ) i n R- 

All tensor products in this work will be over R and all ^-modules 
(resp. ideals) will be left modules (resp. left ideals). 

2.2. Local cohomology. It turns out that Hj(M) may be computed 
as follows. Let C'(fi) be the complex — ► R — * Rf t — > 0. Then 
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Hf(M) is the k-th cohomology group of the Cech complex defined by 
C'{M; f\,... ,f r ) = (g>l C'(fi) ® M. Unfortunately, explicit calcula- 
tions are complicated by the fact that Hf(M) is rarely finitely gener- 
ated as -R-module. This difficulty disappears for Hf(R) if we enlarge 
the ring to A n , in essence because Rf is finitely generated over A n for 
all feR. 

2.3. .D-modules. A good introduction to D-modules is the book by 
Bjork, @. ' 

Let / G R. Then the -R-module Rf has a structure as left A n -module: 
x i (fi) = *ft,d i (f;) = giM^ffiWg . This may be thought of as a special 
case of localizing an A n -module: if M is an A^-module and / G R then 
Rj®rM becomes an y4 n -module via di{j^®m) = di{j^)®m+-j^®dim. 
Localization of A n -modules lies at the heart of our arguments. 

Of particular interest are the holonomic modules which are those 
finitely generated y4 n -modules N for which Ext J An (N, A n ) vanishes un- 
less j = n. Holonomic modules are always cyclic and of finite length 
over A n . Besides that, if N = A n /L, f G R, s is an indeterminate 
and g is some fixed generator of N, then there is a nonzero polynomial 
b(s) in K[s] and an operator P(s) G A n [s] such that P(s)(f ■ f s <S> g) = 
b(s) ■ f s (8> g- The unique monic polynomial that divides all other poly- 
nomials satisfying an identity of this type is called the Bernstein poly- 
nomial of L and / and denoted by bj(s). Any operator P(s) that 
satisfies P(s)f s+1 <S> g = bj(s)f s <8> g we shall call a Bernstein operator 
and refer to the roots of bj(s) as Bernstein roots of / on A n /L. 

Localizations of holonomic modules at a single element (and hence at 
any finite number of elements) of R are holonomic (see ffl, section 5.9) 
and in particular cyclic over A n , generated by f~ a g for sufficiently large 



a G N (see also our proposition [4. 2|) . So the complex C'(N; fi, . . . , f r ) 



consists of holonomic y4 n -modules whenever is holonomic. This facil- 
itates the use of Grobner bases as computational tool for maps between 
holonomic modules and their localizations. As a special case we note 
that localizations of R are holonomic, and hence finite, over A n (since 
R = A n /A is holonomic). 



2.4. The Cech complex. In [11] it is shown that local cohomology 



modules are not only D-modules but in fact holonomic: we know al- 
ready that the modules in the Cech complex are holonomic, it suffices 
to show that the maps are A^-linear. All maps in the Cech complex are 
direct sums of localization maps. Suppose Rf is generated by f s and 
Rf g by (fgY- We may replace s, t by their minimum u and then we see 
that the inclusion Rf Rf g is nothing but the map A n j ann(/ u ) — > 
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A n j ann((/g)") sending the coset of the operator P to the coset of the 
operator P ■ g\ So C\N; f u . . . , f r ) -> C i+1 (N; f u . . . , f r ) is an A n - 
linear map between holonomic modules for every holonomic N. One 
can prove that kernels and cokernels of ^4 n -linear maps between holo- 
nomic modules are holonomic. Holonomicity of Hf(R) follows. 

In the same way it can be seen that H^Hj^R)) is holonomic for 
i,j£N (since Hj(R) is holonomic). 

3. Grobner bases of modules over the Weyl algebra 

In this section we review some of the concepts and results related 
to the Buchberger algorithm in modules over Weyl algebras. It turns 
out that with a little care many of the important constructions from 
the theory of commutative Grobner bases carry over to our case. For 
an introduction into non- commutative monomial orders and related 
topics, || is a good source. 

Let us agree that every time we write an element in A n , we write it 
as a sum of terms c a px a d^ in multi-index notation. That is, a, (3 G W 1 , 
c a p are scalars, x a = x" 1 ■ . . . • x° n , d@ = df 1 ■ . . . • d% n and in every 
monomial we write first the powers of x and then the powers of the 
differentials. Further, if m = c a pX a dP, c a p G K, we will say that m has 
degree degm = |a + (3\ and an operator P G A n has degree equal to 
the largest degree of any monomial occuring in P. 

Recall that a monomial order < in A n is a total order on the mono- 
mials of An, subject to m < m' mm" < m'm" for all nonzero 
monomials m,m',m". Since the product of two monomials in our no- 
tation is not likely to be a monomial (as diXi = Xidi + 1) it is not 
obvious that such orderings exist at all. However, the commutator 
of any two monomials mi,m 2 will be a polynomial of degree at most 
degm! + degm 2 — 2. That means that the degree of an operator and 
its component of maximal degree is independent of the way it is repre- 
sented. Thus we may for example introduce a monomial order on A n 
by taking any monomial order on A n = K[x\, . . . , x n , d\, . . . ,d n ] (the 
polynomial ring in 2n variables) that refines the partial order given by 
total degree, and saying that mi > m 2 in A n if and only if m\ > m 2 in 
A n . 

Let < be a monomial order on A n . Let G = ® x A n ■ ji be the free 
A n -module on the symbols 71, . . . , jd- We define a monomial order on 
G by mi'-fi > mj'jj if either > mj in the order on A n , or mj = mj 
and i > j. The largest monomial rrry in an element g G G will be 
denoted by in(g). Of fundamental importance is 
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Algorithm 3.1 (Remainder). Let h and g = {gi}{ be elements of G. 
Set ho = h, o"o = 0, j = and let e, = e(gi) be symbols. Then 



Repeat 

If in(<7i)| in(hj) set 



Sh — h in ( fc j) n 

_ _i_ wy f 

<7 J+ 1 .- (Tj + 



3 :=J + 1} 
Until Noin(^j)| in(^). 

The result is h a , called a remainder dt(h,g) of h under division by g, 
and an expression o a = Y2i=i a i £ i with a,i e A n . By Dickson's lemma 
([[J, 1.1) the algorithm terminates. It is worth mentioning that dt(h, g) 
is not uniquely determined, it depends on which we pick amongst 
those whose initial term divides the initial term of hj. 

Note that if h a is zero, a a tells us how to write h in terms of g. Such 
a a a is called a standard expression for h with respect to {g\, ... ,g s }. 

Definition 3.2. If in(^) and in(^) involve the same basis element of 
G, then we set = m^gi — m^gj and o~ij = m^Ei — TUijEj where 
rriij = lcm (" 1 i ^' in ( g ')) . Otherwise, and are defined to be zero, s^- 
is the Schreyer-polynomial to $ and grj. 

Suppose ?R(sij,g) is zero for all i, j. Then we call g a Grobner basis 
for the module A n • (pi, ... , g s ). 

The following proposition (|J, Lemma 3.8) indicates the usefulness 
of Grobner bases. 

Proposition 3.3. Let g be a finite set of elements of G. Then g is a 
Grobner basis if and only if h G A n g implies 3i : in(pj)| in(/i). □ 

Computation of Grobner bases over the Weyl algebra works just as 
over polynomial rings: 

Algorithm 3.4 (Buchberger). Input: g = {g±, . . . ,g s } C G. 
Output: a Grobner basis for A n ■ (gi, ... ,g s ). 
Begin. 

Repeat 

If h = ^( Slj ,g) ^0 

add htog 

Until all3?(sy,g) = 0. 
Return g. 

End. 
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3.1. Now we shall describe the construction of kernels of ^-linear 
maps using Grobner bases. Again, this is similar to the commutative 
case and we first consider the case of a map between free y4 n -modules. 

Let E = 0* A n Ei, G = 0^ An^j and : E — > G be a A n -linear map. 
Assume </>(£*) = ft. Suppose that in order to make g a Grobner basis 
we have to add g[, . . . ,g' s , to g which satisfy g\ = J2l=i a ik9k- We get 
an induced map 

S+S ' A p- 

where it is the identity on £j for i < s 

and sends £ i+s into X}fc=i a %k^k- Of course, = 07r. 

The kernel of is just the image of the kernel of under ir. So in 
order to find kernels of maps between free modules one may assume 
that the generators of the source are mapped to a Grobner basis and 
replace by 0. Recall from definition [3.2| that cr^ = m^Ei — rriijSj or 
zero, depending on the leading terms of ft and gj. 

Proposition 3.5. Assume that {gi, . . . ,g s } is a Grobner basis. Let 
Sij = dijkdk be standard expressions for the Schreyer polynomials. 
Then {a^ - J2k dijk£k}i<i<j<s generate the kernel of : 0* A n Ei -> 
0i A n -fj, sending e { to ft. 

The proof proceeds exactly as in the commutative case, see for ex- 
ample ||, section 15.10.8. 




3.2. We explain now how to find a set of generators for the kernel 
of an arbitrary A^-linear map. Let E, G be as in subsection [3.1| and 
suppose A n (pi,... ,p a ) = P C E,A n (q l ,... ,q b ) = Q C G and : 
0i A n 6i/P — > 0iv4„7j/Q. It will be sufficient to consider the case 
P = since we may lift to the free module E surjecting onto EjP. 

Let as before 0(£j) = ft. A kernel element in E is a sum a^j, G 
A n , which if £j is replaced by ft can be written in terms of the generators 
qj of Q. Let j3 = {Pi, . . . , f3 c } be such that gUqU (3 is a Grobner basis 
for A n (g,q). We may assume that the $ are the results of applying 
algorithm |3.4] to gU q. Then from algorithm I3TT] we have expressions 



(3.1) 
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with Cij,c' ik G A n . Furthermore, by proposition p~5| , algorithm |3T4 
returns a generating set g_ for the syzygies on g U q U (3. Write 

Oi = Y a iJ £ 9i + Y a 'ik £ i k + Y a il £ & ( 3 - 2 ) 

3 k I 

and eliminate the last sum using the relations ( |3.1| ) to obtain syzygies 

&i = Yl a ^ £ 9j + "Ys a 'ik £ <lk + Y2 ( ° lv£ 9v + YL C ' lw£ <l™ I • 

j k I \ v w J (3.3) 

These will then form a generating set for the syzygies on gUq. Cutting 
off the g-part of these syzygies we get a set of forms 



Y + Y a ii (Y1 c iv £ 9^j | 



which generate the kernel of the map E —>■ G/Q. 

3.3. The comments in this subsection will find their application in al- 
gorithm ^]2| which computes the structure of H^H^R)) as yl n -module. 
Let 



M' 2 



M 2 - 



MS 



Ml 



M[ —U. Mi — - Ml 

be a commutative diagram of A n -modules. Note that the row coho- 
mology of the column cohomology at N is given by 

ker(V0 n /?' _1 (imp) + im(^)j / [/3(ker(0')) + im(V>)] • 

In order to compute this we need to be able to find: 

• preimages of submodules, 

• kernels of maps, 

• intersections of submodules. 

It is apparent that the second and third calculation is a special case 
of the first: kernels are preimages of zero, intersections are images of 

preimages (if A n r — > A n s /M is given, then im(0) D im('0) = 

V# _1 M$)) ). 

So suppose in the situation : A n r /M — > A n s /N, ip : A^ / P — > 
A n s /N we want to find the preimage under ip of the image of </>. We 
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may reduce to the case where M and P are zero and then lift 0, ip to 
maps into A n s . The elements x in -?/> _1 (im0) C A n l are exactly the 

elements in ker(A„* — > A n s /N — > A^ 8 / (AT + im0)) and this kernel can 
be found according to the comments in |3.2| . 



4. P-MODULES AFTER KASHIWARA AND MALGRANGE 

The purpose of this and the following section is as follows. Given 
/ G R and an ideal L C A n such that A n /L is holonomic and L is 
/-saturated (i.e. / • P G L only if P G I), we want to compute the 
structure of the module Rf <g> A n /L. It turns out that it is useful to 
know the ideal J L {f s ) which consists of the operators P(s) G A n [s] 
that annihilate / s <8> 1 G M := Rf[s]f s Cg> A n /L where the bar denotes 
cosets in A n /L. In order to find J L {f s ), we will consider the module 



M over the ring A n+1 = A n (t, d t ). It will turn out in |4.1| that one can 
easily compute the ideal J J f +1 (/ s ) C A n+1 consisting of all operators 
that kill f s ® 1. In section [5] we will then show how to compute J L (f s ) 
from J% +1 (f s ). 

The second basic fact in this section (proposition ^4.2| ) shows how 
known. 



to compute the structure of Rf £g> A n /L as A„-module once J L (f s ) is 



4.1. Consider A n+1 = A n (t,d t ), the Weyl algebra in x 1; . . . ,x n and 
the new variable t. B. Malgrange has defined an action of t and d t on 
M = R f [s] ■ f s ® R An/L by t{g(x, s) ■ f s ®P)_= g(_x, s + l)f ■ f s ®P 
and d t (g(x, s) ■ f s ® P) = ^#(x, s — 1) ■ f s P for P G A n /L. A n acts 
on M as expected, the variables by multiplication on the left, the d{ 
by the product rule. 

One checks that this actually defines an structure of M as a left 
A n+ i-module and that —dtt acts as multiplication by s. 

We denote by J^ +l {f s ) the ideal in A n+ i that annihilates the element 
f s (g> 1 in M. Then we have an induced morphism of A n+ i-modules 
A/JZ +1 (f s ) - M sending P + J n L +1 (/ s ) to P(f s ® T). 

The operators i and ^ were introduced in Jl3| . The following lemma 
generalizes lemma 4.1 in |13| (as well as part of the proof given there) 
where the special case L = (d%, . . . , d n ), A n / L = R is considered. 

Note that Jn+i(f s ) makes perfect sense even if L is not holonomic. 

Lemma 4.1. Suppose that L = A n ■ (Pi, . . . , P r ) is f -saturated (i.e., 
if f ■ P G L, then P G L). With the above definitions, Jn+i{f s ) is the 
ideal generated by f — t together with the images of the Pj under the 
automorphism <fi of A n+ i induced by x — > x,t —> t — f . 
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Proof. The automorphism sends <9j to <9j + fid t and d t to d t . So if we 
write Pj = Pj(d u . . . , d n ), then <j>Pj = P 5 (di + fid u . .. ,d n + f n d t ). 

One checks that (<9j + fid t )(f s ® Q) = f s <S> 9jQ fc> r all differential 
operators Q so that 4>{Pj{d x , ... , d n )){f s ®T) = f^P^, . . . ,d n ) = 0. 
By definition, / • (f s <g> I) = t ■ (f s ® T). So t - / £ ^+i(/ s ) and 
^)6J n y/ s )for ! =_l 1 ... ) r. 

Conversely let P(f s (g) 1) = 0. We may assume, that P does not con- 
tain any t since we can eliminate t using f—t. Now rewrite P in terms of 
d t and the d { + Say, P = £ c a pd?x p Q a p(di + fidt, . . . ,d n + f n d t ), 
where the Q a p are polynomials in n variables and c a p £ i^. Application 
to / s ® T results in ]T <9f (/ s ® c a px^Q a p(di, ... ,d n )). 

Let a be the largest a £ N for which there is a nonzero c a| g occuring in 
P = S c af3 d?x l3 Q a p{dx + fid t , ... , d n + f n dt). We show that the sum of 
terms that contain df is in A n+ \ -4>{L) as follows. In order for P(/ S (g>i) 
to vanish, the sum of terms with the highest s-power, namely s a , must 
vanish, and so J2p c ap(-^/f) a f s ® x^Qap{di, ... , <9 n ) £ P// s <g> £ as 
Rf is P-flat. It follows, that Yp c a/3 x ^Qoci3(9i, ■ ■ ■ > ^n) £ (-^ is /- 
saturated!) and hence J2p d f c aP^Qap(di+fid t , ■ ■ ■ ,d n +f n d t ) £ A„, +r 

So by the first part, P - Yp Caffi % P Q<xp{<h + fid t , ■■■ , d n + f n d t ) 
kills f s ® 1, but is of smaller degree in <9 4 than P was. 

The claim follows. □ 



4.2. Let J L (/ S ) stand for the ideal in A n [s] = A n [-d t t] that kills 
f s <g> I £ P/[ S ]/ S ® R A n /L. Note that J L (/ S ) = J^ +1 (f s ) n A n [-<9 t t]. 
Again, we may talk about J L (f s ) independently of the holonomicity of 
L. 

We will in the next section show how the lemma can be used to 
determine J L (f s ). Now we show why J L (f s ) is useful, generalizing 



10|| , proposition 6.2. 

Recall that the Bernstein polynomial bh (s) is defined to be the monic 
generator of the ideal of polynomials b(s) £ K[s] for which there exists 
an operator P(s) £ A n [s] such that P(s)(f s+1 <g> T) = b(s)f s <g> T (0, 
chapter 1), and that 6f (s) will exist for example if L is holonomic. 

Proposition 4.2. If L is holonomic and a £ Z is such that no integer 
root ofbj(s) is smaller than a, then we have isomorphisms 

R f <g> AJL = A n [s}/J L (f)\ s=a = A n ■ f a ®l. (4.1) 

Proof. We mimick the proof given by Kashiwara, who proved the propo- 
sition for the case L = (di, . . . , d n ), A n j L = R ( ||10| , proposition 6.2). 
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Let us first prove the second equality. Certainly J L (f s )\ s=a kills 
f a <g> T. So we have to show that if P(f a ® T) = then P G J L {f s ) + 
A n [s] ■ (s — a). To that end note that st acts as t(s — 1) which means 
that t-(A n [s\/J L (f s )) is a left A„[s]-module. Identify A n [s]/J L (f s ) with 
Nf := A„[s]-(/ s ®T). By definition, bj(s) is the minimal polynomial for 
which there is P(s) with ^(s)(/ s ®T) = P(s)f s+1 = t-P(s-l)(/ s ®I). 
So 6y(s) multiplies A n [s] • (/ s <8> 1) into i • A„[s](/ S Cg> 1) and whenever 
the polynomial b(s) G K[s] is relatively prime to bf(s) its action on 
Nf/t ■ Mf is injective. 

Since by hypothesis s — a + j is not a divisor of 6^ (s) for < j G N, 

(s - a + j)Mf nt-Nf Q(s-a + j)t ■ Mj. (4.2) 

So - a + m)JV/ n t m AT f L C (s - a + m)W/ n fW/ = t[(s - a + m - 
1)^/ n t m_1 A/'/'] whenever m > 1. 

We show now by induction on m that (s - a + m)A/} L fl t m Nf C 
(s — a + m)t m Aff for m > 1. The claim is clear for m = 1 from equation 
( |4.2| ). So let m > 1. The inductive hypothesis states that (s — a + m — 
1)^ n t m ~ x N]- C (s — a + m — \)t m ~ x Nj. The previous paragraph 
shows that (a - a + m)A/} L fl fW/ Ct[(«-a + m- 1)-/V/ n t m_ W/] . 
Combining these two facts we get 

{s-a + m)N}nt m Nf C t(s - a + m - l)t m "W/ 

= {s-a + m)t m Nf. 

Now if P(s) G AJs] is of degree m in the 9, and P(a)(f a ® 1) = 0, 
then P(s + m)- f m + J L (f s ) G (s — a+m)-J\ff because we can interprete 
P(s + m)(f s+m (g) 1) as a polynomial in s + m with root a. But then 
P(s + m) (/ s+m (8)1) = P(s + m) (/"7 s <g> T) is in 

(s - a + m)Mf n t m A/} L C (s - a + m)t m Mf, 

implying P(s + m)(f s+m ® 1) = (s — a + m)Q(s)(f s+m ® 1) for some 
Q(s) G A„[s] (note that J L (/ S ) kills / s (g) I). In other words, P(s) - 
(s — a)Q(s — m) G J L (f s )- 

For the first isomorphism we have to show that A n • (/ a <g> 1) = 
Rf <g> A n /L. It suffices to show that every term of the form f m f a <8> <5 
is in the module generated by (f a Cg> 1) for all m G Z. Furthermore, we 
may assume that Q is a monomial in <9i, . . . ,d n . 

Existence and definition of bj(s) provides an operator P(s) with 
[bj(s - l)] _1 P(s - l)(/ s <g> T) = /~7 S ® T. As b L f (a - m) ^ for 
all < m G N we have f m f a <g> T G A n • (/ a <g> T) for all m. Now 
let Q be a monomial in di, . . . , d n of <9-degree j > and assume that 
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f m f a (g> Q' G A n - (f a (g) 1) for all m and all operators Q' of <9-degree 
lower than j. Then Q = diQ' for some 1 < i < n. Fix m G Z. By 
assumption on j, for some P' we have P'(f a <E> 1) = / m / a <8> Q'- So 

f7 a ® Q = ^P'(/ a <g> T) - /i • (a + m)f m - l f a <g> Q 7 g A n • (/« ® T). 

(4.3) 

The claim follows by induction. This completes the proof of the propo- 
sition. □ 

We remark that if any a G Z satisfies the conditions of the proposi- 
tion, then so does every integer smaller than a. 



5. An algorithm of Oaku 

The purpose of this section is to review and generalize an algorithm 
due to Oaku. In |L4| (algorithm 5.4.), Oaku showed how to construct a 
generating set for J L (f s ) in the case where L = (<9i, . . . , d n ). According 
to |3, J L {f s ) is the intersection of J„ +1 {f s ) with A n [-d t t}. We shall 
explain how one may calculate J (1 A n [—d t t] whenever J C A n+ i is any 
given ideal and as a corollary develop an algorithm that for /-saturated 
A n /L computes J L (f s ). The proof follows closely Oaku's argument. 

On A n+1 [y 1 ,y 2 \ define weights w(t) = w(y 1 ) = l,w(d t ) = w(y 2 ) = 
-l,w(xi) = w(di) = 0. If P = Y^iPi e A n+1 [y 1 ,y 2 ] and all P t are 
monomials, then we will write (P) h for the operator ^ Pj ■ yf* where 
di = maxj(w(Pj)) — w(Pi) and call it the yi-homogenization of P. 

Note that the Buchberger algorithm preserves homogeneity in the 
following sense: if a set of generators for an ideal is given and these 
generators are homogeneous with respect to the weights above, then 
any new generator for the ideal constructed with the classical Buch- 
berger algorithm will also be homogeneous. (This is a consequence 
of the facts that the y\ commute with all other variables and that 
dtt = tdt + 1 is homogeneous of weight zero.) 

Proposition 5.1. Let J = A n+ i ■ (Qi, . . . ,Q r ) and let yi,y 2 be two 
new variables. Let I be the left ideal in A n+ i[yi] generated by the y%- 
homogenizations (Qi) h of the Qi, relative to the weight w above, and 
let I = A n+1 [yi,y 2 ] • (1,1 — y\y 2 ). Let G be a Grobner basis for I 
under a monomial order that eliminates yi,y 2 . For each P G G set 
p> = t -w(P)p -y w (pj < q md p, = Q™( p )p ^ w (p) > o and let 

G' = {P' : P EG}. Then G = G' n A n [-d t t) generates J n A n [-d t t). 



Proof. Note first that G consists of w-homogeneous operators and so 
w(P) is well defined for P G G. 
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Suppose P G G . Hence Pel. So P = Q_ x ■ (1 -y x y 2 ) + E a i ■ (<3i) h 
where the a* are in A„ + i[t/i, 7/2]- Since P G A n [— t t], the substitution 
Vl -> 1 shows that P = 1) ■ 1) = E«i(!> !) ' Q< e J. 

Therefore G C J n A n [-d t t}. 

Now assume that P G JnA„[— ftt]. So P is 777-homogeneous of weight 
0. Also, P G J and J is contained in 1(1), the ideal of operators Q(l) C 
A n+ i for which Q(yi) G J. By lemma |5.2| below (taken from 
7/£P G I for some a G N. Therefore P = (1 - (7/i7/ 2 ) a )P + (7/i7/ 2 ) a P G J. 

Let C7 = {Pi, ... , p,, Pft+i, . . . , P c } and assume that p G A n+1 if 
and only if i < b. Buchberger algorithm gives a standard expression 
P = ^2 ciiPi with all in (di Pi) < in(P). That implies that a^ +i is zero 
for positive i and does not contain 7/1,7/2 f° r an Y 

Since P, Pi are 777- homogeneous, the same is true for all ai, from 
Buchberger algorithm. In fact, w(P) = + t/7(P) for all i As 

tt;(P) = (and t, ft are the only variables with nonzero weight that 
may appear in ai) we find a- G A n with = a l i -t~ w< ^ p ^ or a* = a^-d™^, 
depending on whether w(Pi) is negative or positive. 

It follows that P = E? a,P = a[P[ G A„[-ftt] • G . □ 



Lemma 5.2. Lei I be a w -homogeneous ideal in A n+ i[yi] with respect 
to the weights introduced before the proposition and 1(1) defined as in 
the proof of the proposition. Assume P G A n+ i is a w -homogeneous 
operator. Then P G 1(1) implies y®P G / for some a. 

Proof. Note first that 7/1 —>■ 1 will not lead to cancellation of terms in 
any homogeneous operator as 777(7/1) 7^ 0. 

If P G 1(1), P = ^Qi(l), with all Qi 70-homogeneous in I. Then 
the 7/1-homogenization of Qi(l) will be a divisor of Qi and the quotient 
will be some power of 7/1, say y^ . Homogenization of the equation 
P = ^Qi(l) results in y\P = J2Qi(^) h (since P is homogeneous) so 
that 

yV+maxivi) p = J2 y™ x{m) - m Qi G /. □ 

So we have 

Algorithm 5.3. Input: / G R, L C A n such that L is /-satuarated. 
Output: Generators for J L (f s ). 
Begin 

1. For each generator Qi of L compute the image 4>(Qi) under Xi — > 
Xi,t -> t - f,di -> ft + /ift, ft -> ft. Add t - / to the list. 

2. Homogenize all <p(Qi) with respect to the new variable 7/1 relative 



to the weight 77; introduced before proposition O. 
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3. Compute a Grobner basis for the ideal generated by (<j>(Qi)) h , 

(0(Qr))\ 1 - yW2, t - yif in A n+1 [yx, y 2 ] using an order that 
eliminates y\,yi- 

4. Select the operators {Pj}\ in this basis which do not contain yi, y 2 . 

5. For each Pj, 1 < j < b, if w(P j ) > replace Pj by Pj = dt (Pj) P r 
Otherwise replace Pj by Pj = t~ w ( p ^Pj. 

6. Return the new operators {Pj}i- 
End. 

In order to guarantee existence of the Bernstein polynomial b p (s) we 
assume for our next result that L is holonomic. 

Corollary 5.4. Suppose L is a holonomic ideal. If J L {f s ) is known or 
it is known that L is f -saturated, then the Bernstein polynomial b^(s) 
ofR f ® R A n /L can be found from (tf(s)) = A n [s] ■ (J L (/ S ), /) H K[s). 

Moreover, if K C C, suppose bj(s) = s d + bd~is d ~ l + . . . + b and 
define B = maxj{|6j| 1 ^ d_ *- ) }. In order to find the smallest integer root 
ofbj(s), one only needs to check the integers between —2B and 2B. 

If in particular L = (di, . . . ,d n ), it suffices to check the integers 
between —bd~i and -1. 



Proof. If L is /-saturated, propositions [O] and |5.1| enable us to find 
J L (f s ). The first part follows then easily from the definition of bj(s): 
as (bj(s) - Pj ■ f)(f s ® T) = it is clear that b L f {s) is in K[s\ and 
in A n [s](J L (f s ), f). Using an elimination order on A n [s], bj(s) will be 
(up to a scalar factor) the unique element in the reduced Grobner basis 
for J L (f s ) + (f) that contains no 2j nor 

Now suppose K C C, |s| = 2Bp where B is as defined above and 
p > 1. Assume also that s is a root of bj(s). We find 

d-l d-1 

(2Bp) d =\s\ d = | -J2 biSl \ ^ J2 Bd ^\ s \ l t 5 ' 1 ) 



d-l 



B d J2(2pY<B d ((2p) d -l), (5.2) 



using p > 1. By contradiction, s is not a root. 

The final claim is a consequence of Kashiwara's work |TI| where it is 
proved that if L = (d\, . . . , d n ) then all roots of bj(s) are negative and 
hence — & n -i is a lower bound for each single root. □ 

For purposes of reference we also list algorithms that compute the 
Bernstein polynomial to a holonomic module and the localization of a 
holonomic module. 
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Algorithm 5.5. Input: / 6 R, L C A^ such that A n /L is holonomic 
and /-torsionfree. 

Output: The Bernstein polynomial bj(s). 

Begin 

1. Determine J L (f s ) following algorithm B^B. 



2. Find a reduced Grobner basis for the ideal J L (f s ) +A n [s] ■ f using 
an elimination order for x and d. 

3. Pick the unique element in that basis contained in K[s] and return 
it. 

End. 



Algorithm 5.6. Input: / G R, L C A n such that A n /L is holonomic 
and /-torsionfree. 

Output: Generators for an ideal J such that Rf £g> A n /L = A n /J. 

Begin 

1. Determine J L (f s ) following algorithm |5.3| . 

2. Find the Bernstein polynomial bj(s) using algorithm |5.5| . 

3. Find the smallest integer root a of bj(s) (using corollary if 
KCC). 

4. Replace s by a in all generators for J L (f s ) and return these gen- 
erators. 

End. 



The algorithms [5^3] and [5.5| appear already in |14[ in the special case 
L — (<9i, . . . , d n ),A n / L — R. 



6. Local cohomology as A^-module 

In this section we will combine the results from the previous sections 
to obtain algorithms that compute for given i,j,k&N,lQR the local 
cohomology modules Hj(R), H^Hji^R)) and the invariants Xij(R/I) 
associated to /. 

6.1. Computation of Hj(R). Here we will describe an algorithm that 
takes in a finite set of polynomials / = {fi, . . . , f r } C R and returns a 
presentation of Hj(R) where I = (/i, . . . , f r ). In particular, if Hf(R) 
is zero, then the algorithm will return the zero presentation. 
Consider the Cech complex associated to fi, . . . , f r in R, 

r 

<' •0/'/ ' "/J %..,/. "<>• 

1 l<i<j<r (6-1) 
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Its k-th cohomology group is the local cohomology module H k (R). The 
map 

C k = Rf h .....f ih -> R fji ..... fjk+i = C k+1 

l<h<-<i k <r l<ji<-<jk+i<r (6.2) 

is the sum of maps 

R f . . f . -> i? f . . . f . (6.3) 

which are either zero (if . . . ,4} {ji> • • ■ jjfc+i}) or sen d j to 
y, up to sign. Recall that A n /A = A n /A n ■ (di, . . . ,d n ) = R and 
identify />'./>..,./,. with A n /J A ((f h f ik ) s )\ s =a and /.'/,,.,/„. with 
An/ J ((f fa • . . . ■ /j fe+1 ) s )| s =f, where a, 6 are sufficiently small integers. 
By proposition |4.2| we may assume that a = b < 0. Then the map (|6.2[) 



is in the nonzero case multiplication from the right by (/;) a where 
l = {ji,... ,jk+i}\ {h,--- again up to sign. 

It follows that the matrix representing the map C k — > (7 fc+1 in terms 
of A n -modules is very easy to write down once the annihilator ideals 
and Bernstein polynomials for all k- and (k + l)-fold products of the 
fi are known: the entries are or ±ff a where fi is the new factor. 

Let Q r k be the set of fc-element subsets of 1, . . . , r and for 8 G Q r k write 
F g for the product riiee r /»■ We have demonstrated the correctness and 
finiteness of the following algorithm. 

Algorithm 6.1. Input: /i, . . . , f r e R; k 6 N. 

Output: Hj(R) in terms of generators and relations as finitely gen- 
erated y4„-module. 
Begin 

1. Compute the annihilator ideal J a ((Fq) s ) and the Bernstein poly- 
nomial bp g (s) for all (k — 1)-, k- and {k + l)-fold products of 
/i s , as in U and |J (so 6 runs through Q r k _ x U0 r fc U Q r k+1 ). 

2. Compute the smallest integer root ag for each bp* (s), let a be the 
minimum and replace s by a in all the annihilator ideals. 

3. Compute the two matrices M k _i,M k representing the y4 n -linear 
maps C k ~ l — > C fc and C — > as explained in subsection |Q . 

4. Compute a Grobner basis G for the kernel of the map 

©4-©4/J A TOU- A„/j A ((F e )% =a 



as in |3T2| . 

5. Compute a Grobner basis Go for the module 

im(Mfc_i) + J A ((F e ) s )\ s=a C A n /J A ((F 
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6. Compute the remainders of all elements of G with respect to lifts 
of C7 to 0ee r A n . 

7. Return these remainders and Gq. 

End. 

The nonzero elements of G generate the quotient G/G = Hj(R) so 
that Hj(R) = if and only if all returned remainders are zero. 



6.2. Computation of H l m {H 3 j{R)). As a second application of Grob- 
ner basis computations over the Weyl algebra we show now how to 
compute H^Hj^R)). Note that we cannot apply lemma pL"T| to A n /L = 
Hj(R) since H^R) may well contain some torsion. 

As in the previous sections, C^(R; fi, . . . , f r ) denotes the j-th mod- 
ule in the Cech complex to R and {fi, . . . , f r }. Let C" be the double 
complex with C lJ = C l (R; x±, . . . , x n ) <S>r C 3 (R; fi, . . . , f r ), the verti- 
cal maps 0" induced by the identity on the first factor and the usual 
Cech maps on the second, whereas the horizontal maps £** are in- 
duced by the Cech maps on the first factor and the identity on the 
second. Since C l (R; x±, . . . ,x n ) is i?-projective, the column cohomo- 
logy of C" at is C l (R; Xi, . . . ,x n ) Hj(R) and the induced 

horizontal maps in the j-th row are simply the maps in the Cech com- 
plex C m (Hj(R); xi, . . . , x n ). It follows that the row cohomology of the 
column cohomology at (io,Jo) is H^°(Hf(R)). 

Now note that C M is a direct sum of modules R g where g = x ai ■ 
x ai ■ fp 1 ■ . . . ■ fp, . So the whole double complex can be rewritten 
in terms of A„-modules and yl n -linear maps using |5\6| : 

Qi-lj+V >. (Ji,j+1 — — >. 



i— 1,3 



»-l,J-l 



Using the comments in subsection |3.3| , we may now compute the mod- 
ules H^HjlR)). More concisely, we have the following 

Algorithm 6.2. Input: fi, . . . , f r G R; io,jo G N. 

Output: H^(Hf (R)) in terms of generators and relations as finitely 
generated A n -module. 

Begin. 
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1. For i = i — l,i ,i + 1 and j = j — l,j ,j + 1 compute the 
annihilators J A ((Fg ■ Xg/) s ) and Bernstein polynomials bp g , x t (s) 
of Fg ■ Xgi where 6 e 6' G 0™ and Xg> denotes in analogy to 
F e the product Ylaee' x »- 

2. Let a be the minimum integer root of the product of all these 
Bernstein polynomials and replace s by a in all the annihilators 
computed in the previous step. 

3. Compute the matrices to the j4 n -linear maps JJ : C l ' J — > C % ' 3+ 
and : G M — > C* +1,J ', again for z = i — 1, z , «o + 1 an d j = 

Jo ~ 1, jo, jo + 1. 

4. Compute Grobner bases for the modules 

G = kei{<f) io ' jo ) n [(f °' jo )" 1 (im(^ 0+lj °- 1 ))] + im(0 ioJ '°- 1 ) 

and G = f o " 1 ' jo (ker(0 io - ljo )) + im(^°' J '°- 1 ) . 

5. Compute the remainders of all elements of G with respect to Go 
and return these remainders together with Go- 

End. 

The elements of G will be generators for H^(Hj°(R)) and the ele- 
ments of Go generate the relations that are not syzygies. 

6.3. Computation of Xi tn -j(R/I). In [I1J h has been shown that 



H^HjlR)) is an injective m-torsion i?-module of finite socle dimen- 
sion \i, n -j (which depends only on i,j and R/I) and so isomorphic 
to (En(K)) Xi < n -i where E R (K) is the injective hull of K over R. We 
are now in a position that allows computation of these invariants of 
R/I. For, let Hl l (H 3 I (R)) be generated by g\, . . . , gi £ A n d modulo 
the relations h\,... ,h e £ A n d . Let H be the module generated by the 
hi. We know that (A n • g\ + H) /H is m-torsion and so it is possible 
(with trial and error) to find a multiple of g\, say mgi with m a mono- 
mial in R, such that (A n ■ mg\ + H ) / H is nonzero but Xitngi e H for 
all 1 < i < n. Then the element mgi + H/H has annihilator equal 
to m and hence generates an ^-module isomorphic to A n /A n ■ m = 
E R {K). The injection A n ■ mg x + H/H A n • . . . + H/H 
splits as map of i?-modules because of injectivity and so the cokernel 
A n (g 1 , . . . ,gi)/An(mgi, h, . . . , h e ) is isomorphic to {E R {K)) Xi > n -i~ l . 

Reduction of the gi with respect to a Grobner basis of the new rela- 
tion module and repetition of the previous will lead to the determina- 
tion of \i, n -j. 

6.4. Local cohomology in ambient spaces different from A^. 

If A equals K[xi, . . . ,x n ], I C A, X — Spec(A) and V = Spec(A/I), 
knowledge of H}(A) for all i e N answers of course the question about 
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the local cohomo logical dimension of V in X. It is worth mentioning, 
that if W C X is a smooth variety containing V then our algorithm 



6.1 for the computation of H}(A) also leads to a determination of the 
local cohomological dimension of V in W. Namely, if J stands for the 
defining ideal of W in X so that R = A/ J is the affine coordinate ring 
of W and if we set c = ht(J), then it can be shown that H]~ C (R) = 
Rom A (R, Hj(A)) for all i G N. As Hj(A) is /-torsion (and hence J- 
torsion), Hom^(i2, Hj(A)) is zero if and only if H}(A) = 0. It follows 
that the local cohomological dimension of V in W equals cd(A, I) — c 
and {q G N : Hj(A) ^ 0} = {q G N : Hj~ c (R) ^ 0}. 

If however W is not smooth, no algorithms for the computation of 
either H){R) or cd(R,I) are known, irrespective of the characteristic 
of the base field. 

7. Implementation and examples 

Some of the algorithms described above have been implemented as 
C-scripts and tested on some examples. 

7.1. The algorithm [5.3| with L = A has been implemented by Oaku 



using the package Kan (see JL7]]) which is a postscript language for 
computations in the Weyl algebra and in polynomial rings. An imple- 
mentation for general L is written by the current author and part of a 
program that deals exclusively with computations around local coho- 
mo logy (fTJfl). [0 is theoretically able to compute H}(R) for arbitrary 
i,R = Q[xi, . . . , x n ], I C R in the above described terms of generators 



and relations, using algorithm |6J]. It is expected that in the near future 
18 1 will work for R = K[xi, . . . , x n ] where K is an arbitrary field of 
characteristic zero and also algorithms for computation of H^Hj^R)) 
and Xij(R) will be implemented, but see the comments in [7.2| below. 



Example 7.1. Let / be the ideal in R = K[xi, ... ,xq] that is gener- 
ated by the 2x2 minors f , q, h of the matrix I Xl X2 X3 ] . Then 

\ x 4 x 5 x e J 

Hj(R) = for i < 2 and Hj(R) ^ because I is a height 2 prime 
and H}(R) = for i > 3 because I is 3-generated, so the only remain- 
ing case is Hf(R). This module in fact does not vanish, but until the 
discovery of our algorithm, its non-vanishing was a rather non-trivial 
fact. Our algorithm provides the first known proof of this fact by direct 
calculation. 

Previously, Hochster pointed out that Hf(R) is nonzero, using the 
fact that the map K[f,g,h] — > R splits (compare 0, Remark 3.13) 
and Bruns and Schwanzl (0, the corollary to Lemma 2) provided a 
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topological proof of the nonvanishing of Hf(R) via etale cohomology. 
Both of these proofs are quite ingenious and work only in very special 
situations. 

Using the program fL8f , one finds that Hf(R) is isomorphic to a cyclic 



Ag-module generated by 1 6 A 6 subject to relations x 1 — . . . — x e — 0. 
This is a straightforward computational proof of the non-vanishing of 
Hf(R). Of course this proof gives more than simply the non- vanishing. 
Since the quotient of A§ by the left ideal generated by 
known to be isomorphic as an i?-module to Er(R/(xi, . . . , x$)), the 
injective hull of Rj (x±, . . . , xq) = K in the category of .R-modules, our 
proof implies that Hf{R) = E R {K). 

7.2. Computation of Grobner bases in many variables is in general a 
time- and space consuming enterprise. Already in (commutative) poly- 
nomial rings the worst case performance for the number of elements in 
reduced Grobner bases is doubly exponential in the number of variables 
and the degrees of the generators. In the (relatively small) example 
above R is of dimension 6, so that the intermediate ring A n+ x[yi,y 2 ] 
contains 16 variables. In view of these facts the following idea has 
proved useful. 



The general context in which lemma [4.1| and proposition |47| were 
stated allows successive localization of Rf g in the following way. First 
one computes Rf according to algorithm |5^ as quotient of A n by a 



certain holonomic ideal L = J A (f s )\ s=a , a <C 0. Then Rf g may be 
computed using |5]6| again since Rf g = R g ® A n /L. (Note that all lo- 
calizations of R are automatically /-torsion free for / e R as R is a 
domain.) This process may be iterated for products with any finite 
number of factors. Note also that the exponents for the various factors 
might be different. This requires some care as the following situations 
illustrate. Assume first that —1 is the smallest integer root of the Bern- 
stein polynomials of / and g (both in R) with respect to the holonomic 
module R. Assume further that Rf g = A n ■ f^g^ 1 2 A n ■ (fg)^ 1 . 
Then Rf —>■ Rf g can be written as A n j ann(/ _1 ) — > A n j ann(/~ 2 ■ g^ 1 ) 
sending P 6 A n to P ■ f ■ g. 

Suppose on the other hand that we are interested in H](R) where 
I — (/; 9-, h) an d we know that Rf = A n - f~ 2 D A n ■ R g = A n ■ g~ 2 
and Rfg = A n ■ f^g^ 2 . (In fact, the degree 2 part of the Cech complex 
of example |7.1| consists of such localizations.) It is tempting to write 
the embedding Rf — > Rf g with the use of a Bernstein operator (if 
Pf(s)f s+1 = bf~(s)f s then take s = —2) but as f^ 1 is not a generator 
for Rf, bf L (—2) will be zero. In other words, we must write Rf g as 
A^/ ann((/g)™ 2 ) and then send P G ann(/~ 2 ) to P ■ g 2 . 
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The two examples indicate how to write the Cech complex in terms 
of generators and relations over A n while making sure that the maps 
C k — ► C k+1 can be made explicit using the ff the exponents used in 
C % have to be at least as big as those in C 1 " 1 (for the same /j). 

Remark 7.2. We suspect that for all holonomic /g-torsionfree mod- 
ules M = A n /L we have (with R g ® M = A n /V): 

min{s G Z : bj(s) = 0} < min{s G Z : bf (s) = 0}. 

This would have two nice consequences. 

First of all, it would guarantee, that successive localization at the 
factors of a product does not lead to matrices in the Cech complex 
with entries of higher degree than localization at the product at once. 
Secondly, if were known to be true, we could proceed as follows for 



the computation of C l (R; fx, . . . , / r ). First compute J A {{fi) s ) for all 
i, find all minimal integer Bernstein roots $ of /$ on R and substitute 
them into the appropriate annihilator ideals. If from now on we want to 
use algorithm |5.6| in order to compute Rf il -...-f ik -f ik : from Rfi^...-f ik then 
we can skip steps 2 and 3 of as the remark gives us a lower bound 
for the minimal integer Bernstein root of f% k+1 on Rf il ....-f i • (From the 



comments before [Z.2| it is also clear that we cannot hope to use a larger 
value.) 

The advantage of localizing Rf g as (Rf) g is twofold. For one, it allows 
the exponents of the various factors to be distinct which is useful for 
the subsequent cohomology computation: it helps to keep the degrees 
of the maps small. (So for example R x .( x 2 +y 2^ can be written as A n ■ 
x^ 1 (x 2 + y 2 )~ 2 instead of A n ■ (x~ 2 ■ (x 2 + y 2 Y 2 )- On the other hand, 
since the computation of Grobner bases is doubly exponential it seems 
to be advantageous to break a big problem (localization at a product) 
into many "easy" problems (successive localization). 

An extreme case of this behaviour is our example \! .1\ if we compute 
Rfgh as ((Rf) g )h, the calculation uses approximately 2.5 kB and lasts 
32 seconds on a Sun workstation using [[[J|. If one tries to localize R 
at the product of the three generators at once, [jn| crashes after about 
30 hours having used up the entire available memory (1.2 GB). 
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